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I.  CEPSTRUM  ANALYSIS 

(J.  A.  Shooter,  C.  J.  Webb) 

It  was  pointed  out  in  the  previous  progress  report^  that  the 
techniques  for  computing  a Cepstrum  have  certain  shortcomings  in  a 
mathematical  sense.  These  restrictions  were  placed  on  the  technique 
by  certain  approximations  and  assumptions  that  were  made  concerning 
the  relative  amplitudes  of  the  original  signal  and  echoes.  This 
past  quarter,  effort  was  made  to  derive  alternate  techniques  for 
computing  Cepstrum,  using  as  few  mathematical  assumptions  as  possible 


In  the  following  discussions,  the  standard  method  for  computing 

2 

a Cepstrum  as  defined  by  B.  P.  Bogert  et  al,  will  be  referred  to  as 
Cepstrum.  The  alternate  methods  will  be  defined  as  derived. 

A.  First  Alternate  for  Cepstrum  Analysis  (DRL  A- I 


Mathematical  Description  for  a Single  Echo 


Consider  the  time  series 


z(t)  = y(t ) + Ay(t  - t)  , 


Defense  Research  Laboratory,  The  University  of  Texas  at  Austin, 
Quarterly  Progress  Report  No.  3 under  Naval  Ship  Systems  Command 
Contract  NObsr-95l8l,  1 February  - 30  April  1967,  pp.  2-7. 

B.  P.  Bogert,  M.  J.  R.  Healy,  and  J.  W.  Tukey,  in  "Proceedings 
of  the  Symposium  on  Time  Series  Analysis,"  ed.  M.  Rosenblatt 
(John  Wiley  and  Sons,  Inc.,  New  York,  1963),  Chap.  15,  pp.  209-243. 
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consisting  of  the  signal  y(t)  together  with  the  same  signal  delayed 
in  time  by  the  amount  t and  damped  in  amplitude  by  the  factor  A. 


The  following  notation  will  be  used  throughout: 

fz(^)  = Fourier  transform  of  z(t), 
fy(tu)  = Fourier  transform  of  y(t), 

$ (u>)  = power  spectrum  of  z(t), 

F (o>)  = power  spectrum  of  y(t). 


o>  = angular  frequency  in  rad/sec. 


Then,  from  Eq.  (l). 


*(“>  ■ f>>  • ■ fy(“)  • fy%)(l  + Aeta)(l  * Ae-toT)  , (; 

where  the  asterisk  denotes  complex  conjugation.  Therefore, 

<5(03)  = F(cu)(l  + 2A  cos  cut  + A2)  . 

For  those  values  of  cu  for  which  *(<u)  / 0,  F(cu)  / 0,  *»  and  F'(cu) 
exist ; 


m F 'U} 
H<x>)  F(cu ) 


t sin  cut 


(A  + l)/2A|  f cos 


Thus,  the  logarithmic  derivative  of  the  power  spectrum  of  the 
composite  signal,  z(t),  is  the  logarithmic  derivative  of  the  original 
signal,  y(t),  with  an  added  periodic  function. 


variable?  den°teS  * deriVative  wlth  to  the  dependent 


- t sin  cot 


(A"  + 1)/2A 


+ COS  COT 


This  added  tern  is  periodic  in  co/2n  with  period  1/t . Hence,  if 


CO 


is  relatively  flat,  this  additional  term  should  appear  as  a 


"ripple"  of  period  1/t  superimposed  on  the  log  derivative  of  F(co). 


$ ^ ((i)  ) 

If  a spectral  analysis  of  is  now  performed,  and  if 

is  relatively  flat,  the  component  with  frequency  1/t  should 
contribute  significantly.  Specifically,  examining  the  power 

. . (J)  f ( /J\  \ 

spectrum,  P(a),  of  a peak  at  a = t is  expected. 


In  contrast  to  Cepstrum,  this  approach,  to  be  referred  to  as 
DRL  A- I,  requires  no  linear  approximation  to  a logarithm  function 
Also,  no  restrictions  have  been  placed  on  the  relative  amplitudes 
of  the  delayed  and  undelayed  signals  as  is  necessary,  at  least 
theoretically,  in  the  case  of  Cepstrum. 


2.  Practical  Comparison  of  Cepstrum  and  DRL  A- I 

The  alternate  method  for  computing  a Cepstrum,  DRL  A-I, 
is  not  more  difficult  to  implement  on  the  digital  computer  than  the 
standard  Cepstrum  technique.  The  logarithmic  derivative  can  be 
computed  by  a standard  difference  scheme. 

As  an  example  of  the  methods,  consider  the  signal  z(t)  to 
contain  only  one  echo  as  described  by  the  relation 

z(t)  = sin[2rtft]  + H(t  - t)  sin[2rtf(t  - t)]  , (5) 
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given  as 


(6) 

t > 1 . 

The  sample  length  will  be  2.0  sec  with  a bandwidth  from  50  to  150  Hz. 
The  sampling  rate  is  1000  Hz. 

Figure  1 shows  the  time  function  given  by  Eq.  (5)  processed  by 
both  the  Cepstrum  and  DRL  A- I techniques.  As  expected,  both  methods 
give  the  same  results  for  this  simple  model. 

3.  Multiple  Echoes 

Consider  first  the  case  of  two  time  delays  or  echoes 

given  as 


where  the  amplitude  of  the  echo  is 


z(t)  = Y(t)  + A1y(t  - tx)  + A2y(t  - T ) . (7) 


In  this  case, 


f,(cn)  = fy(u>) 


-law.  -iu>T 

1 + A1e  + A^ 


, and  (8) 


k 


®((o)  = F(co) 


- IcOT  - IcOT 

1 + A^e  + Ag° 


/ 2 p 

= F(co)  S 1 + + A^  + 2A1  cos  cot 


(9) 


Thus, 


+ 2a^  cos  cot ^ + 2A^2  COS  a)(Ti  “ T2M 


j>/(co)  _ F/(co) 

4>(co)  _ F(<o) 

A1t1  sin  <ot  + A2Tg  sin  coTg  + A^  - Tg)  sin  o>(t  - t ) 
5 J5  “ 

1 + A1+  A2 


(10) 


2 + A±  COS  COT1  + A2  COS  COT2  + A1A2  cos  co(t2  - T ) 

^ an<^  ^2  are  sma-^  m°<lulus  compared  to  unity,  and  the  higher 

order  terms  in  A^  are  negligible,  then 


tt'fcp)  = F'(co)  A1T1  stn  “Ti  + A2t2  Sin  m2 
®(“)  F(co)  l/2  + A cos  cot  + A cos  cot,. 

12  c. 


(n) 


The  numerator  and  denominator  of  the  second  term  on  the  right  side 
of  Eq.  (11)  each  consist  of  terms  periodic  in  co/2t  with  periods 
l/T1  and  1/t2. 

In  the  case  of  N delays,  it  is  simple  to  see  that 


6 


’",,p  " ’V'"'T'A 


*'{*)  3 F/(aO 

^(dt)  f(o>) 


1/2  + 


£ 


A,  cos  ayr, 
k k 


(12) 


It  is  of  interest  to  note  that  the  simplification  recommended 
(ignoring  the  higher  order  terms  in  the  A's)  can  be  made  either 
before  or  after  taking  the  logarithmic  derivative;  i.e.,  either  in 
El • (9)  °r  in  Eq.  (10).  Equation  (ll)  is  the  result  regardless. 

4.  Single  Echo  With  Noise 

If  we  consider  a background  noise  in  addition  to  the 
signal,  the  simple  model  described  by  Eq.  (l)  becomes 

z(t)  = n(t ) + y(t)  + Ay(t  - t)  . (13) 

In  this  case 

*z(»)  = Fn(a>)  + 2Re  jfjoo)  f *(o>)  (l  + Ae^) 

(14) 

+ E^(u>)  [l  + A"1  + 2A  cos  uuj  , 

where  the  asterisk  indicates  complex  conjugation,  Re  denotes  the 
real  component  of  the  expression  in  brackets,  and  the  other  notation 
is: 

f (u>)  = Fourier  transform  of  n(t), 
fy(co)  = Fourier  transform  of  y(t), 


l 

) 
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Fn(a>)  = Power  spectrum  of  n(t), 
F^(a>)  = Power  spectrum  of  y(t). 


Equation  (14)  shows  the  degree  to  which  noise  can  compound  the 
problem.  However,  if  n(t)  is  a white  noise  source — clearly  the  most 
desirable  situation — then  F (u>^  is  flat.  F /(u))  therefore  vanishes 
and  the  Fn(u>)  term  makes  no  contribution  to  the  logarithmic  deriva- 
tive of  $(o)).  Note  that  the  second  term  on  the  right  hand  side  of 
Eq.  (l4)  will  tend  to  zero  with  increasing  values  of  signal  to  noise 
ratio.  Thus,  in  the  case  of  a white  noise  source  and  high  signal - 
to-noise  ratio,  the  logarithmic  derivative  of  4>z(u>)  [4>z(u.')  as  given 
by  Eq.  (l4)]  can  be  made  to  approximate  Eq.  (4). 


B.  Second  Alternate  for  Cepstrum  Analysis  (DHL  A- II 


From  the  previous  development,  it  is  clear  that  a signal, 
y(t),  with  an  optimally  flat  power  spectrum  is  the  most  desirable 
for  additional  processing  such  as  Cepstrum  or  DRL  A- I.  For  example, 
if  the  signal  is  a sine  pulse,  the  frequency  is  modulated  to  achieve 
a power  spectrum  that  is  approximately  flat . 


A flat  power  spectrum  can  be  obtained  by  a more  judicious 
choice  of  the  signal  y(t).  Suppose  that  y(t)  is  a signal  whose 
Fourier  transform,  f^(<n),  possesses  the  properties: 


f (u>) 

yv  ’ 


B,  0 < u>  < cn 

* - - o 


0,  to  < to 
o 


and  f (u>)  is  real. 


For  a signal  composed  of  y(t)  and  a single  echo,  such  as 


z(t)  = y(t)  + y(t  - t)  , 


the  Fourier  transform  is 


f»  =< 


B j_l  + Ae"  , 0 < cd  < 


CU 


0 


,0)  <01 

o 


(16) 


(IT) 


The  Imaginary  component  of  fz(oo)  is 

Im  fz(^)  = -AB  sin  cut  . (l8) 

Thus,  Im  fz(u>)  consists  of  a single  periodic  component,  and  the  power 
spectrum  of  Im  f“z (oo ) is  a delta  function  of  time. 

For  the  case  of  multiple  echoes  or  N time  delays,  it  can  be 
shown  that  the  result  will  be  N delta  functions- -one  delta  function 
at  each  of  the  points  in  time  t^,  k = 1,  2,  N. 

The  advantages  of  this  second  alternate  method  lie  in  its 
obvious  simplicity  and  short  computation  time.  However,  whereas 
Cepstrum  and  DRL  A- I are  designed  for  use  with  a more  general  type 
signal,  this  last  method  requires  a signal  with  very  special 
properties.  One  example  of  a signal  of  this  type  is 

y(t)  - 2^  , (i9) 

which  has  a rectangular  shaped  Fourier  transform. 
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II.  DIGITAL  FILTERS 

(G.  E.  Ellis,  H.  W.  Pitman) 

A.  Nonrecursive  Filters 

For  a linear  filter  described  by  the  weighting  function 
h(t),  the  output  of  the  system  u(t),  can  be  related  to  the  system 
input,  u(t),  over  finite  time  limits  TL,Tu  by  the  relation 


For  discrete  or  sampled  data,  this  expression  can  be  given  as 

u 

g = y h u , (21) 

D?  c_j  n n4  in 
n=-L 

where  h are  the  filter  weights,  u , are  the  input  samples,  and 
n n+m 

u are  output  samples . The  spacings  n and  m are  determined  by  the 
m 

sampling  rate  needed  to  prevent  aliasing. 

For  infinite  time  limits  (-«  < n < °°)  the  frequency  function, 
H(o>),  of  the  filter  together  with  the  weighting  function,  continuous 
or  discrete,  form  particular  Fourier  transform  pairs . In  the  case 
of  finite  time  limits,  it  is  desirable  to  choose  a finite  set  of 
weights,  hn,  whose  frequency  is  a least  square  fit  to  a designed 
frequency  function. 

A class  of  filters  having  a frequency  response  function  of 
unity  gain  and  zero  phase  shift  in  the  passband  with  a pth  order 
roll-off  between  the  passband,  (to  < <x>c ) , and  the  rejection  bonds, 

(cjd  is  given  in  Fig.  2. 
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FIGURE  2 

FREQUENCY  FUNCTION  FOR  A CLASS  OF  LINEAR  FILTERS 

0»c  = 27Tfc  = CUTOFF  FREQUENCY 

t Oj  : 2-rrfj  = FILTER  ROLL-OFF  TERMINATION  FREQUENCY 


DRL  - UT 
AS. 68-43 
GEE  - EJW 
1 - 26  - 68 


H(a>)  can  be  specified  as 


H(a>)  = ( 


0 

1 


00  > 03l 


M < % 


("T  ' -c)‘  ^ + ‘V  ' “ ^T  - 


< d)  < -03 


^ (aVp  - cc)-P  (a)  - a^,)p,  a>c  < a»  < a>T 


(22) 


In  the  time  domain,  the  weighting  function  is  given  by 


(23) 


In  terms  of  the  number  of  weights  for  a numerical  filter,  a higher 
order  of  decay  implies  less  weights  for  a given  accuracy  in  the 
desired  frequency  representation  of  the  finite  filter.  The  expres- 
sion for  the  weights  increases  rapidly  with  increasing  order  of 
decay.  If  use  is  made  of  the  fast  Fourier  transform  (FFT)  for  the 
convolution,  Eq.  (21 ),  large  numbers  of  weights  are  practical  and 
extremely  sharp  filter  responses  can  be  attained. 

For  the  case  of  p = 1,  the  weighting  function  is  given  as 

COS  0)  t - COd  Gli_t 

h(t)  - — -£ -P-  • (24) 

The  discrete  weights,  h , can  be  evaluated  directly  from  h(tn)At, 

*\j.  F.  A.  Ormsby,  "Design  of  Numerical  Filters  with  Applications  to 
Missile  Data  Processing,"  J.  ACM,  8,  440-466  (1961). 
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where  equidistant  samples  are  assumed.  The  assumption  that  h(t)  is 
real  gives  H(-cd)  = H*(cd),  where  * denotes  the  complex  conjugate. 

For  zero  phase  shift,  the  frequency  response  function  is  even  or 
H(u>)  = |h(oo)  | and  h = h_n-  In  Eq.  (2l),  u = L = N,  to  give  a total 
of  2N  + 1 weights. 

For  ease  of  analysis,  the  frequency  will  be  normalized  by 
the  sampling  frequency  cDg  to  give  a nondimens ional  variable  r = CD/cDg, 
Equation  (24)  then  becomes 


cos(2nnr  ) - cos(2nnrrT,) 
h_  = 


2r  (*n)£ 


n = 0,  ±1,  ...,  ±N  , (25) 


where 


rc  ■ “A 
rT  " V“S 

CD  = 2nf 


r = r - r 
r T c 

t = nAt  , 
n ’ 


h = h(t^)At  , and 


At  = 2jt/cd 


The  relation  for  the  zero  order  weight  is 


h = r_  + r . (26) 

o T c 


These  relations  need  to  be  subjected  to  certain  constraints^  for 

d£ 


CD. 


-*■  0,  especially  that  H(cd)  = 1 and 


12H^cd)  = 
dcD 


0 at  cd  = 0.  To 


Ibid. 
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determine  the  constrained  set  of  filter  weights,  the  following 
relations  are  used: 


h 

o 


= h + 
o 


ctA  + 2 cl  5, 


2 1 


12 


(27) 


h = h + g26l  + 2glS2  ~ ^ + (2N  + 1)8J 

n n a " > n > 1 , (28) 

where 


n=l 


a5  = (2N  + l)o2  - 2a,2 


For  a high-pass  filter  with  similar  characteristics,  the 
weights  (hQ)  are  determined  frcm  the  low- pass  weights  (h  ) by  the 
relations  n 
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.Jjr  7W3T-A 


h = 1.0  - h 
o o 


h = -h 
n n 


(29) 


A band-pass  filter  can  be  derived  by  subtracting  the  weights  of  two 
low-pass  filters,  but  for  symmetry  a preferable  method  is  frequency 
shifting.  If  the  low- pass  filter  frequency  function  is  given  by  H(r) 
and  weights  hQ,  the  frequency  function  for  the  band-pass  filter  is 
defined  as  Hg(r)  = H(r  - rQ)  + H(r  + rQ),  where  rQ  is  the  center  of 
the  desired  bandpass.  For  zero  phase  shift  filters. 


r* 

yr)  - 2\  ♦ 2 £ 
n=l 


(2E  cos  2jtnrQ)  cos  2jtnr 


(30) 


The  weights  for  Kg(r)  are  given  by 


hn  = 2hn  COS  (2ltnr0) 


(31) 


Quadrature  filters  or  filters  with  90  deg  phase  shift  in 
the  passband  can  be  derived  in  a similar  manner.  Once  again,  it  is 
more  expedient  to  derive  a low-pass  filter  first  and  then  convert  to 
a band-pass  or  high-pass  filter.  Accuracies  corresponding  to  the 
zero  phase  shift  filters  can  be  obtained. 


Figure  3 shows  the  relation  between  accuracy  over  the  range 
0 - f ^ fc>  the  number  of  filter  weights,  and  rp,  a parameter  that 
controls  the  sharpness  of  the  roll-off.  Figure  4 gives  the  relation 
between  sampling  rate,  the  number  of  filter  weights,  and  the  band- 
width of  a band-pass  filter.  The  inherent  problem  with  nonrecursive 
filters  is  apparent  from  these  curves--to  construct  a narrow-band 
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FIGURE  4 

NUMBER  OF  WEIGHTS  = (2N  + 1) 

DRL  - UT 
AS.  68-45 

BW  = FILTER  BANDWIDTH 

GEE  - EJW 
1 - 26  - 68 

filter  at  sampling  rates  near  10  KHz,  it  is  necessary  to  have  a large 
number  of  filter  weights.  This  problem  can  be  offset  scmewhat  by 
using  the  fast  Fourier  transform  algorithm  for  performing  the  convo- 
lution with  the  data  series. 

B.  Recursive  Filters 


Recursive  filtering,  or  more  simply  a filter  that  operates 
on  a data  series  with  a recursive  relationship,  is  analagous  to  an 
analog  filter  with  feedback.  The  digital  application  of  filters  of 
this  type  is  enhanced  since  a recursive  filter  normally  is  derived 
with  only  a few  weights.  The  convolution  or  filtering  process  does 
not  demand  large  blocks  of  computing  time  as  in  the  case  of  non- 
recursive filters. 

An  example  of  the  design  of  recursive  filters  with  the  use 
of  Z-transformations  will  now  be  outlined.  The  filter  to  be  con- 
sidered is  a high-pass  RC  filter.  The  circuit  and  gain  character- 
istics are  shown  in  Fig.  5. 

The  first  step  in  the  digital  simulation  is  to  describe 
the  continuous  transfer  function  of  the  filter  by  Laplace  transform 
techniques.  Since  samples  from  the  continuous  data  are  assumed  c,o 
be  taken  at  discrete  intervals  of  time  a,  a unit  delay  is  repre- 

a 

sented  by  multiplying  the  continuous  Laplace  function  by  e , 
where  s is  the  continuous  Laplace  variable.  The  unit  delay  interval 
can  be  expressed  in  the  manner  of  the  Z-transform  as 


Z = e~as  = 1 - as  + (as)2/2'.  - (as)^/3*  + ...  , (32) 

or 

s = -|(logeZ)  . (33) 
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a 


It  is  time-consuming  to  substitute  for  s in  this  form.  By  taking  an 
approximation  to  Z,  the  algebra  is  simplified  and  the  results  still 
match  the  analog  filter  results  very  closely.  The  approximation  for 


Z = 1 - as  + (as)2/2  - (as)^/2  • 2 + (as)^/2  • 2 ' 2 - ...  , (34) 


which  reduces  to 


- 1^1 

a U * S| 


The  steps  for  converting  an  analog  filter  to  a digital  recursive 
filter  by  this  technique  can  be  summarized  as: 

l)  Put  the  analog  filter  transfer  function  in  Laplace  transform 


form; 

2)  Replace  the  Laplace  s by  WN 
sampling  time  of  the  data  series; 


where  At  is  the 


3)  Reduce  the  resulting  equation  to  the  form 


u 

-Xa 


B ZJ  E . 
j out 


Interpret  Z in  Eq.  (36)  to  be  a unit  sample  time  delay 
(Z  = two  sample  times,  etc.); 

4)  Write  the  equation  for  the  digital  filter  in  the  form 


i 


[•  - 


•>  1 wmomm 


m 


"out 


<*>  ' I A1  Eln  <r  - l>  - I B<  E„„t  <*  - J> 


(37) 


i=o 


J=1 


For  the  RC  high- pass  filter  the  steps  are  as  follows: 


(1) 


out 


RCs 


E.  1 + RCs 
in 


(38) 


(2) 


out 


M 

L1  * 


in 


1 + A 


1_  7 
~ 

1 + z 


(39) 


where  A = 2RC/At; 


(3)  E = E [(1  - Z)  t-j-t]  - E [z  ] 

out  in  ' 1 + AJ  out  (_  1 + A J ; 


m 


[ 


W EOUt  ^ = n (I)  - E<n  (I  - 1) 


in 


in 


)1  \ - A 

J l 1 + A 


- Eout 


i + aJ 


W 


Figure  6 shows  the  frequency  response  of  the  RC  high-pass 
filter  given  by  Eq.  (4l).  As  in  the  analog  case,  the  ideal  curve 
cannot  be  attained.  The  digital  filter  simulates  the  analog  filter 
extremely  well,  with  a -3  dB  point  at  l/RC  or  0.1  Hz,  and  a slope 
of  -6  dB/octave.  If  the  user  wishes  a filter  with  steeper  slopes, 
the  filter  may  be  used  repeatedly  or  in  cascade.  An  example  of 
cascading  is  given  in  Fig.  7*  The  filter  was  used  twice,  or  in  the 
analog  sense  as  a two- stage  RC  filter,  to  give  a -6  dB  point  at 
l/RC  and  a -12  dB/octave  slope. 
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IDEAL  FILTER 
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IDEAL  FILTER 


FIGURE  7 

RESPONSE  OF  CASCADED  RC  HIGH-PASS  FILTER 


In  terms  of  computing  time,  a filter  of  this  type  is 
extremely  efficient  and  fast.  It  can  be  cascaded  many  times  for 
very  steep  slope  without  a great  expenditure  of  computer  time . The 
digital  simulation  of  other  analog  filters  can  be  accomplished  in  a 
similar  manner.  Several  very  good  references  on  the  subject  of 
recursive  digital  filters  are: 

1)  Blackman,  R.  B.  and  Tukey,  J.  W.  The  Measurement  of  Power 
Spectra  (Dover  Publication,  New  York,  i960) . 

2)  Golden,  R.  M.  and  Kaiser,  J.  F.  "Design  of  Wideband 
Sampled- Data  Filters,"  Bell  System  Tech.  43,  1533-1646, 

(1964). 

3)  Jury,  E.  I.  Theory  and  Application  of  the  Z- Trans form 
Method  (John  Wiley  and  Sons,  Inc.,  New  York,  1964). 

4)  Kuo,  F.  F.  and  Kaiser,  J.  F.  System  Analysis  by  Digital 
Computer  (John  Wiley  and  Sons,  Inc.,  New  York,  1966), 

218-277. 

5)  Shanks,  John  L.  "Recursion  Filters  for  Digital  Processing," 
Paper  presented  at  the  36th  Meeting  of  the  Society  of 
Exploration  Geophysicists,  Houston,  Texas,  1967  (Unpublished). 
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III.  DIGITAL  SIMULATION  LANGUAGES 
(G.  E.  Ellis) 


The  discrete  simulation  language  GASP  II  has  been  implemented 
on  the  CDC  3200  computer.  GASP  (General  Activity  Simulation  Program) 
consists  of  thirteen  subroutines  and  functions  written  in  FORTRAN 
and  specifically  designed  for  use  in  simulation  studies  of  industrial 
systems. 

GASP  II  provides  the  programmer  with  a library  of  subprograms 
which  may  be  used  to  accomplish  certain  operations  common  to  many 
simulation  studies  and  models.  In  general,  GASP  will  perform  the 
following  functions:  start  the  simulation,  move  the  model  through 
time,  stop  the  simulation,  monitor  and  present  intermediate  results 
of  the  simulation  if  desired,  provide  random  variables  drawn  from 
the  more  common  probability  distributions  according  to  parameters 
specified  by  the  programmer,  collect  and  classify  data,  compute 
statistics,  detect  and  isolate  errors  in  the  program,  and  present  a 
summary  report  of  the  simulation. 

The  present  version  of  GASP  II  was  written  for  a computer  with 
a limited  version  of  FORTRAN  and  a small  memory.  This  quarter  has 
been  devoted  to  updating  the  coding  with  the  more  flexible  FORTRAN 
available  for  the  CDC  3200  and  expanding  the  program  size  to  fit  the 
32k  memory  bank  of  the  3200. 
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IV. 


PERT 

TlTa.  Noack,  G.  E.  Ellis) 

A model  research  program  was  devised  to  include  experimental, 
theoretical,  and  data  reduction  phases.  The  program,  to  be  con- 
ducted during  a one- year  period,  is  shown  in  Fig.  8.  The  network 
is  arranged  in  a normal  PERT/TIME  configuration.  The  critical  path, 
or  longest  path  in  time,  is  shown  as  derived  from  the  computer 
output  from  the  CDC  3200. 

At  first  glance,  the  outstanding  disadvantage  that  appears 
in  the  network  is  the  lack  of  a feedback  loop,  which  is  essential 
in  all  research  programs.  Using  the  standard  PERT /TIME  computer 
programs,  the  feedback  loop  can  be  implemented  by  the  scientist 
and  the  computer,  updating  with  new  projected  times  and/or  a modi- 
fication of  the  network. 

The  simulation  language  GASP  II  will  be  used  in  the  next 
quarter  to  implement  PERT/TIME  techniques.  If  this  is  successful, 
a feedback  loop  in  the  process  will  be  implemented  using  GASP  II 
to  more  closely  model  the  actual  research  process. 


29 


F 


f 

■ 


i 


Copy  No. 
1 - 5 


8 

9 

10 

11 

12 

13 

14 

15 

16 

17 


2 February  1968 


DISTRIBUTION  LIST  FOR 
QUARTERLY  PROGRESS  REPORT  NO.  4 
ON  CONTRACT  NObsr-95l8l 
PROJECT  SERIAL  NO.  SF1010316,  TASK  864l 
1 May  - 31  July  1967 


Commander,  Naval  Ship  Systems  Command  (SHIPS  00V1C) 
Department  of  the  Navy 
Washington,  D.  C.  20360 


Director, 

U.  S.  Naval  Personnel  Research  Activity 

Department  of  the  Navy 

San  Diego,  California  92152 


Office  of  Naval  Research 
Resident  Representative 
The  University  of  Texas  at  Austin 
2507  Main  Building 
Austin,  Texas  78712 


Signal  Physics  Division,  DRL/UT 

K.  J.  Diercks,  DRL/UT 
G.  E.  Ellis,  DRL/UT 

R.  K.  Goodnow,  DRL/UT 

L.  A.  Jeffress,  DRL/UT 

S.  P.  Pitt,  DRL/UT 

J.  A.  Shooter,  DRL/UT 
R.  H.  Wallace,  DRL/UT 
C.  J.  Webb,  DRL/UT 
Library,  DRL/UT 


31 


DEFENSE 

RESEARCH 

LABORATORY 


